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SUMMARY 


The  corrections  for  lift,  drag,  and  moment  of  a  two-dimensional  air 
foil  have  been  analyzed,  on  the  assumption  that  the  airfoil  is  tested  in 
the  working  section  of  a  supersonic  tunnel  in  which  the  pressure  field, 
instead  of  being  uniform,  is  characterized  by  gradients  in  the  axial  and 
transverse  directions. 

The  pressure  gradients  of  the  tunnel  as  well  as  the  effect  'of  the 
airfoil  to  be  tested  are  regarded  as  perturbations  of  the  original  rec¬ 
tilinear  flow  field  of  given  Mach  number.  Therefore  the  velocity  poten¬ 
tial  of  the  flow,  the  nonlinear  differential  equation  of  motion,  and  the 
boundary  conditions  are  expanded  into  double  series  in  powers  of  two 
parameters,  one  characterizing  the  airfoil  thickness  e  and  the  other, 
the  inhomogeneity  of  the  field  b.  In  this  way  the  nonlinear  problem 
is  split  into  a  system  of  linear  boundary- value  problems,  corresponding 
to  the  different  powers  of  b  and  e . 

In  each  of  the  resulting  problems  there  appears,  besides  the  dif¬ 
ferential  equation  and  boundary  condition,  an  additional  condition  to 
be  stipulated  on  the  characteristics  passing  through  the  leading  edge. 
Particular  attention  has  been  paid  to  the  correct  formulation  of  this 
"characteristic  condition." 

The  solution  procedure  is  carried  out  up  to  orders  b2 f  e2,  and 

eb  in  the  velocity  potential.  This  means  that,  for  example,  the  drag 

is  computed  up  to  orders  b2e,  and  e2b.  The  physlcal  meaning  of 

the  results  is  discussed.  The  drag  term  in  eb  represents  the 

"horizontal  buoyancy"  of  the  airfoil,  the  term  proportional  to  be2  is 
a  consequence  of  the  interaction  of  the  airfoil  field  and  the  inhomo¬ 
geneous  pressure  field,  and  the  term  in  b2e  may  be  considered  as  a 
"second-order  buoyancy."  The  meaning  of  the  various  lift  and  moment 
terms  may  be  interpreted  similarly.  The  resulting  expressions  have  been 
derived  for  arbitrary  given  pressure  gradients  and  general  profile  form. 

All  solutions  are  obtained  in  closed,  analytic  form  ready  for  imme¬ 
diate  evaluation.  Representative  examples  with  graphs  are  included. 
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INTRODUCTION 


Reports  from  supersonic  tunnels  indicate  that  sometimes  there  exist 
in  the  working  section  undesired  static  pressure  gradients  due  to  design 
or  construction  flaws ,  which  are  difficult  to  eliminate.  As  a  consequence, 
force  measurements  made  upon  models  have  to  be  corrected  to  agree  with  the 
forces  to  be  expected  in  rectilinear  flow. 

For  incompressible  flow  the  correction  to  be  applied  to  the  drag, 
in  the  case  of  a  longitudinal  pressure  gradient  is  well-known  (refer¬ 
ence  l).  Because  of  the  occurrence  of  a  "horizontal  buoyancy,"  the  cor¬ 
rection  turns  out  to  be  equal  to  the  product  of  the  pressure  gradient 
and  the  sum  of  volume  and  "apparent  volume"  of  the  test  body;  the  appar¬ 
ent  volume,  like  the  real  volume,  depends  only  on  the  geometrical  shape. 

The  present  analysis  is  an  attempt  to  solve  the  analogous  problem 
in  supersonic  flow  under  the  most  general  conditions,  which  require  cor¬ 
rections  to  be  applied  also  to  the  lift  and  moment.  A  two-dimensional 
pressure  gradient  in  both  the  longitudinal  a.nd  transverse  directions  is 
considered,  which  is  equivalent  to  a  stream-angle  variation  along  the 
tunnel  axis.  In  the  first  approximation  the  computations  must  be 
expected  to  yield  a  superposition  of  well-known  results.  Hence,  it  is 
important  to  carry  the  calculations  up  to  the  second  approximation, 
which  includes  a  term  characterizing  the  interaction  between  the  airfoil 
field  and  that  of  the  given  pressure  gradient. 

This  work  was  conducted  at  New  York  University  under  the  sponsorship 
and  with  the  financial  assistance  of  the  National  Advisory  Committee  for 
Aeronautics . 


b 

c 

C0 

k(x) 

E(x) 


SYMBOIB  ’ 

parameter  characterizing  inhomogeneity  of  field 
local  sound  speed 
sound  speed  at  Mq 

dimensionless  airfoil  profile  function  measured 
above  flight  direction 

profile  function  measured  above  chord 


Z 


chord  length 
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M 


local  Mach  number 


Mq  Mach  number  of  uniform  field 

m(x),n(x)  x-  and  y-components  of  velocity,  on  axes,  due  to 
inhomogeneity  of  flow  field  (dimensionless) 


P 


local  pressure 


Po 


pressure  at  Mq 


q  local  velocity 

x,y  dimensionless  coordinates,  axial  and  vertical 


3  =  VMo2  -  1 


7 

e 

6 

e 

P 

po 

9 

9* 

9' 

t(x) 


adiabatic  exponent 
airfoil  thickness  parameter 
local  angle  of  attack 
angle  of  attack  of  chord 
characteristic  coordinates 
local  density 
density  at  MQ 

dimensionless  velocity  potential 
velocity  potential  of  basic  flow  field 
velocity  potential  due  to  airfoil  disturbance 
ordinate  of  leading-edge  characteristic 


STATEMENT  OF  PROBLEM 
Physical  Formulation 


Suppose  that  in  the  test  section  of  a  supersonic  tunnel  the  velocity 
field,  instead  of  being  uniform,  is  represented  by  the  velocity  distribu¬ 
tion  along  the  axis  y  =  0,  as  follows: 
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At  y  =  0: 


Mx  =  Mq  +  bm(x) 
My  =  bn(x) 


(1) 


Here  Mx  =  qx/c0  and  My  =  qy/cQ,  where  qx  and  qy  are  the  x-  and 
y-components  of  the  local  velocity  of  the  flow  field,  cQ  is  the  sound 
velocity  in  a  region  where  the  Mach  number  is  Mq,  x  and  y  are 
dimensionless  coordinates  measured  in  units  of  the  length  l  of  the  test 
body,  and  b  is  a  dimensionless  parameter  characterizing  the  magnitude 
of  the  inhomogeneity  of  the  flow  field  and  the  pressure  gradient.  Under 
actual  conditions  b  is  a  small  quantity  such  that 


bm(x) 

bn(x) 


«  Mo 


(2) 


where  m(x)  and  n(x)  are  given  arbitrary  functions  along  the  axis, 

which  may  be  determined  by  measurement.  With  conditions  on  the  axis 

fixed  the  flow  field  in  the  whole  xy-plane  is  automatically  prescribed 

by  the  equations  of  supersonic  flow.  Terms  of  the  order  b2  may  also 

be  given  arbitrarily  on  the  axis  but  are  assumed  to  be  zero  in  the  • 

present  case.  This  does  not  imply,  however,  that  terms  of  the  order  b2 

do  not  appear  in  the  field  away  from  the  axis. 

In  this  velocity  field  a  thin  airfoil  of  given  profile  form  k(x) 
and  thickness  ratio  e  is  inserted,  which  may  be  regarded  as  another 
small  perturbation  of  the  rectilinear  flow  field  Mq.  Two  kinds  of  dis¬ 
turbances  then  exist,  the  one  produced  by  the  velocity  or  pressure  gra¬ 
dient  measured  by  b  and  another  one  produced  by  the  airfoil  and 
measured  by  e.  Therefore,  it  appears  appropriate  to  expand  the  veloc¬ 
ity  potential  cp,  as  well  as  the  potential  equation  of  supersonic  flow, 
and  the  boundary  condition  on  the  airfoil  surface  into  a  power  series  in 
powers  of  both  parameters  b  and  e . 

In  this  way,  the  nonlinear  problem  is  split  into  a  system  of  linear 
boundary- value  problems  corresponding  to  the  different  powers  of  b 
and  e.  In  each  of  the  resulting  problems  there  appears,  besides  the 
differential  equation  and  boundary  condition,  an  additional  condition  to 
be  stipulated  on  the  characteristics  passing  through  the  leading  edge: 

It  is  to  be  required  that  the  disturbance  produced  by  the  airfoil  vanish 
along  these  curved  lines  in  order  to  make  the  solution  unique.  It  may 
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be  noted  that  the  characteristics  of  the  complete  differential  equation 
are  curved,  while  those  of  the  resulting  linear  equations  are  straight. 
It  is  shown  later  how  the  condition  on  the  characteristics  can  be  for¬ 
mulated  for  each  linear  boundary- value  problem  separately. 


Analytic  Formulation 

Because  of  the  relatively  small  thickness  of  the  airfoil,  the  shocks 
at  the  leading  edge  will  be  weak  and  the  entropy  changes  in  the  flow 
negligible.  Hence,  the  complete  flow  field  will  be  described  by  a  veloc¬ 
ity  potential  which  obeys  the  well-known,  nonlinear  differential  equation. 
Introducing  a  nondimens ional  quantity  cp,  equal  to  the  potential  divided 
by  ZcQ,  this  equation  may  be  written: 


{^2  ~  (Px2)(PxX  -  +  ~  V^Cpyy  =  0  (3) 

where  the  subscripts  x  and  y  indicate  differentiation. 

The  local  sound  speed  c  corresponding  to  the  local  velocity  q 
is  related  to  the  sound  speed  cQ  by  means  of  the  energy  equation 


where 

q2/co2  =  'Px2  +  ?y2 


Substituting  equation  (4),  equation 


(3)  can  be  written: 


7+1 

2 


29^9^  + 


(M, 


7+1 


2 


2 


0 


(5) 
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Now,  cp  may  be  expanded  as  a  series  in  powers  of  the  perturbation 
parameters  b  and  e,  as  follows: 

cp  =  cp00  +  bcp01  +  b2cp02  +  bcp10  +  bV1  +  ecp12  +  ebcp13  +  e2cplif  (6) 


The  potential  coefficients 


9®  =  900  +  bcp01  +  b2cp02  (7) 

represent  the  potential  of  the  basic  inhomogeneous  velocity  field  before 
an  airfoil  was  placed  in  it. 

The  terms 


9'  =  b910  +  b291:L  +  e2912  +  eb913  +  g2q)l4 


(8) 


represent  the  perturbation  potential  produced  by  the  airfoil  in  the 
basic  field. 

It  is  to  be  noted  that  the  perturbation  potential  contains  the 

terms  b910  +  b29i;L  which  do  not  vanish  for  an  airfoil  of  zero  thick¬ 
ness  (  e  =  0) . 

A  disturbance  for  which  €  =  0  may  be  thought  of  as  a  very  thin 
flat  plate  placed  at  a  given  angle  into  the  basic  flow  field.  In  order 
for  the  inhomogeneous  flow  to  follow  the  surface  of  the  plate,  vertical 
disturbance  velocities  must  be  produced  of  such  magnitude  (up  or  down¬ 
ward)  as  to  cancel  the  vertical  velocity  components  of  the  inhomogeneous 
flow  field,  which  are  proportional  to  b  or  b2. 

The  potential  term  proportional  to  eb  is  due  to  the  interaction 
between  the  basic  field  and  that  produced  by  the  airfoil.  The  e  and 

e  terms  represent  the  first-  and  second-order  potentials  of  an  airfoil 
in  a  homogeneous  field. 

Finally,  the  analytic  procedure  to  be  followed  in  this  investigation, 
will  be  outlined  briefly: 

(l)  Differential  equation:  Substituting  series  (6)  for  9  into  the 
nonlinear  equation  (equation  (5))  and  splitting  in  powers  of  b  and  e, 
a  set  of  homogeneous  or  nonhomogeneous  equations  of  the  wave  type  is 
obtained. 
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p 

Those  equations  which  are  proportional  to  b,  e,  and  e  are 
homogeneous : 

(»o2  -  iK/*  -  tPyj115  "  0  (9) 

p 

The  equations  which  are  proportional  to  b  and  eb  are  of  the 
nonhomogeneous  type: 


(Mo2  -  l)^xlk  -  ^lk  -  (10) 

where  the  functions  Fjv  can  be  shown  to  be  known  functions  of  lower- 

order  cp^m  which  have  been  computed  before. 

(2)  Boundary  condition:  Insert  series  (6)  for  9  into  the  usual 
boundary  condition  stating  that  the  flow  is  everywhere  tangential  to  the 
surface  of  the  body: 


JPx /surface 


(11) 


and  expand  equation  (ll)  about  y  =  0  in  powers  of  6.  This  imposes 
on  each  of  the  potential  coefficients  cpik  a  condition  of  the  form 


(12) 


where,  again,  the  functions  G-y^  can  be  shown  to  be  known  functions  of 
lower-order  <p^m. 

(3)  "Characteristic  condition":  As  mentioned  before,  the  poten¬ 
tial  9'  describing  the  various  disturbances  produced  by  the  airfoil 
is  required  to  vanish  on  the  leading-edge  characteristics  which  may  be 
written: 


y  =  i(x,b)  =  \lr0(x)  +  bi|/1(x)  +  .  .  . 


(13) 


since  the  form  of  the  characteristics  is  influenced  by  the  inhomogeneity  b. 
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It  can  be  shown  (see  section  "Characteristic  Condition")  that  the 
"characteristic  condition"  for  cp*  is  equivalent  to  the  following 

relations  to  be  imposed  on  the  individual  cp  : 

2 

For  the  orders  b,  e,  and  e 


cplk(x,t0)  =  0 


(ib) 


For  the  orders 


p 

b  and 


be 


cplk(x,\|r0)  =  t1(x)q)^m(x,\(f0) 


(15) 


where,  again,  the  coefficients  cp^m  are  known  lower-order  potentials 
and  i|/q(x)  and  ^(x)  may  be  computed  by  means  of  the  theory  of 

characteristics . 

At  the  end  of  this  procedure  the  solutions  cp^k  will  have  been 
derived  in  closed,  analytic  form,  and  the  pressures  at  any  point  can  be 
computed.  Since  the  objective  is  to  compute  the  corrections  for  forces 
and  moments,  it  is  sufficient  to  determine  the  velocities  and  pressures 
on  the  airfoil  surface  only,  thus  simplifying  the  considerations. 


BASIC  VELOCITY  FIELD 
Differential  Equation 


The  analysis  can  be  clarified  by  considering  first 
geneous  field  alone.  Then  the  series  expansion  for  cp 


the  basic  inhomo - 
reduces  to: 


cp6  =  cp00  +  bcp01  +  b2q)02 


(16) 


2C 
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Substituting  equation  (l)  into  equation  (5)  yields: 


|l  -  Mo2  -  b(r  +  D^cp/1  +  .  .  ^(bcpxx01  +  b2cpXx02  +...)- 


2(Mc  +  .  .  O^bqpy01  +  .  .  •J^xy01  +...)  + 

[l  -  b(7  -  l^qPx01  +  .  .  r|(bcpyy01  +  b2cpyy02  +...)=  0  (17) 


In  deriving  equation  (17)  9^  -  MQx  has  been  substituted  in  equa¬ 

tion  (5),  since  it  is  obviously  a  solution;  furthermore,  all  those  terms 
have  been  kept  which  could  possibly  contribute  to  the  differential  equa¬ 
tion  in  b  and  b2.  Splitting  up  in  powers  of  b,  the  following  equa¬ 
tions  result: 

In  b: 

(m  2  -  l)q>  01  -  9  01  =  0 
\  0  /  hxx  ryy 

P 

and  mb: 

{%2  -  ^'Pxx02  ■  'Pyy02  =  -Mo[(7  +  1^cPx°ltPxx01  +  "  ^Tx^Tyy01  + 

29  019  oil 

yy  yxy  J 

Boundary  Condition 

i 

The  velocity  distribution  (equation  (l)),  prescribed  on  the  axis 
y  =  0,  may  be  considered  as  a  boundary  condition  for  the  potential  of 
the  basic  flow. 
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Substituting  equation  (l6)  into  equation  (l),  it  follows  that 


Mo  +  bcpx01(x,0)  +  b^cp^XjO)  =  Mq  +  bm(x) 


2*  02/ 


and 


b<p  °l(x,0)  +  b2cp  02(x,0)  =  bn(x) 

*/  o 


(18) 


Ordering  terms  in  powers  of  b,  a  boundary' problem  in  b  and  one  in  b2 
are  obtained. 


Boundary- Value  Problem  in  Order  b 

The  boundary-value  problem  in  order  b  consists  of  the  differential 
equation 


p2cp  01  -  cp  01  =  0 
H  'xx  \yy 


(19a) 


O  p 

where  p^  =  MqC  -  1,  and  the  boundary  conditions,  at  y  =  0, 


CPX°1(X^  °)  =  m(x) 


cpy01(x,0)  =  n(x) 


(19b) 


The  solution  of  differential  equation  (19a)  is 


901(x?y)  =  f(x  +  Py)  +  g(x  -  Py) 


(20) 


where  f  and  g  are  arbitrary  functions  to  be  determined  by  boundary 
conditions  (19b). 

Differentiating  equation  (20)  with  regard  to  x  and  y  and  then 
putting  y  =  0  yield: 


01 


cp  (x,0)  =  f' (x)  +  g'(x)  =  m(x) 


cpy01(x,0)  =  p[f'(x)  -  g'  (x)]  =  n(x) 


(21) 


where  f’  and  g'  indicate  derivatives  with  regard  to  the  whole  argument. 
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From  equations  (21)  there  follow: 

f'  (x)  =  ||m(x)  +  i  n(xjj 

g'  (x)  =  |  m(x)  -  -  n(x)J 


(22) 


Since  it  can  he  seen  from  equation  (20)  that  f  and  g  or  their 
derivatives  are  constant  along  their  characteristics  x  ±  Py,  it  may  he 
concluded  from  equations  (22)  that 


and 


f  (x  +  Py)  =  |jm(x  +  fry)  +  ^  n(x  +  py^J 
g' (x  -  py)  =  ||m(x  -  3y)  -  ^  n(x  -  0y)j 


(23) 


Differentiating  equation  (20)  and  substituting  equations  (23),  the 
velocity  components  of  order  b1  for  the  entire  field  follow: 


9x01(x>y)  =  ||m(x  +  py)  +  i  n(x  +  Py)  +  m(x  -  Py)  -  i  n(x  -  pyTJ 
q>y01(x^y)  =  |[m(x  +  0y)  +  |  n(x  +  py)  -  m(x  -  py)  +  ~  n(x  -  PyTJ 


(24) 


From  equations  (24)  the  expressions  for  cpxx<’)1  and  <pXy01  which  appear 
in  the  differential  equations  for  qp-l  and  cp!3  may  he  readily  derived. 

Since,  later  on,  the  values  of  the  velocity  and  pressure  field  are 
required  only  on  the  surface  of  the  body,  that  is,  near  the  axis,  expres¬ 
sions  (24)  may  he  expanded  about  the  axis: 


01(x,ek(x)]]  =  m(x)  +  ek(x)n’ (x)  .  .  . 


cpy01  [x, ek(x)]  =  n(x)  +  ek(x)p2m’  (x)  .  .  . 


(25) 
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Boundary- Value  Problem  in  Order  b2 

In  order  b2  the  boundary-value  problem  consists  of  the  differ¬ 
ential  equation 


and  the  boundary  conditions 

02 

<Px  (x,0)  =  0 
cpy02(x,0)  =  0 


(26a) 


(26b) 


In  order  to  determine  the  velocity  and  pressure  field  on  the  airfoil  sur- 
face ,  no  explicit  solution  of  equations  (26)  is  necessary.  This  can  be 
seen  by  the  following  reasoning. 


For  y  =  ek(x)  the  derivatives  of  <Py02(x,y) 
expanding  about  the  axis: 


may  be  obtained  by 


cpx02|x,ek(x)]  =  cpx02(x,0)  +  €k(x)q>xy02(x,0)  +  .  . 
cpy02£,6k(xri  =  <Py02(x,0)  +  ek(x)cpyy02(x,0)  +  .  . 


V  (27) 


Since  according  to  equation  (l6)  cp02  is  of  order  b2,  the  second-  and 
higher-order  terms  in  equations  (27)  are  of  order  b2e,  b2e2,  .  #  # 
and  may  thus  be  neglected  so  that  in  order  b2: 


<px02  (x,ek(xH  £  0 

Ty02[*,ek(x)]  ~  0 


(28) 
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Therefore,  the  total  velocity  components  of  the  basic  field  on  the 
airfoil  surface  are: 


cpxB  =  Mq  +  bm(x)  +  ebk(x)n’ (x)  +  b20  .  .  . 


cpyB  =  bn(x)  +  ebk(x)p2m' (x)  +  b20  .  .  . 


(29) 


It  is  to  be  noted  that  the  second  derivatives  cpxx®2  and  cpyyP2 

appear  in  the  differential  equation  for  cp11  (see  equation  (32c))  and 
hence  it  may  seem  that  the  boundary- value  problem  (equations  (26))  has 

to  be  solved  after  all.  However,  the  combination  <Pyy02  -  02<PXXO2> 
which  is  required  in  equation  (32c),  may  be  expressed  in  terms  of  the 
known  lower-order  potential  cp01,  as  differential  equation  (26a)  indi¬ 
cates.  Carrying  out  the  proper  differentiations  and  substitutions 
yields: 


where  £  =  x  +  py  and  q  =  x  -  Py  are  the  characteristic  coordinates 
of  the  problem. 
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VELOCITY  FIELD  OF  AIRFOIL 
Differential  Equation 


The  next  objective  is  to  determine  the  perturbation  potential  pro¬ 
duced  by  an  airfoil  which  is  placed  into  the  (inhomogeneous)  basic  field 
at  an  arbitrary  sma.ll  angle  of  atta.ck,  the  leading  edge  coinciding  with 
the  coordinate  origin.  Inserting  series  (6)  for  cp  into  equation  (5), 
the  quantities  cp01  and  cp02  may  now  be  considered  as  known  and 
cp10  •  •  •  Cp^\  as  looked  for.  Keeping  only  those  terms  which  may  pos¬ 
sibly  contribute  to  a  differential  equation  of  the  order  b,  e,  b2. 
be,  or  c  ,  (equation  5)  yields: 


Ordering  terms  in  powers  of  b  and  e,  equation  (31)  will  split  into 
the  following  set  of  linear  equations: 


In  b: 


p 

0  cp. 


10 


xx 


-  cp 


yy 


=  o 


(32a) 
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In  e : 

V  "  ^yy12  =  0  (32b) 

In  b2: 

P^xx11  -  'Pjy11  =  -*kjj5  t  1)  +  (r  -  Dp^'Pxx01  +  ^%x01  *  <px10)  + 
2(£py°1  +  'py10)(9xy°1  +  'Pxy10)}  +  %y°2  -  P2<ft<x02  (32c) 

In  be : 

P29XX13  -  CPyy13  =  -Mq-  [[7  ^  1)  +  (7  -  1)3^1  |^X12(CPXX01  +  CPXXl0)  + 

^xx12^01  +  <Px10j]  +  2^y12(q)xy01  +  cpxy10)  + 

v12!^01  +  V0)]}  (32a> 

P2CPXx14  -  9yy14  =  -Mq  |Tr  +  1)CPX12CPXX12  +  ~  D^Sy^  + 

2V2v12]  (32e) 

Whenever  the  differential  equation  is  inhomogeneous,  the  right-hand  side 
turns  out  to  be  a  function  of  the  lower -order  solutions  which  have  been 
computed  in  an  earlier  step. 

Boundary  Condition 

The  usual  boundary  condition  that  the  velocity  of  the  air  is  every¬ 
where  tangential  to  the  airfoil  surface  may  be  written  as  follows: 
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At  y  =  ek(x) 


Ty  =  ek*(x)cpx  (33) 

As  in  the  basic  field,  so  also  here  the  values  of  the  various  potent ia.l 

coefficients  cp^k  on  the  airfoil  surface  may  be  obtained  by  expansion 
about  y  =  0: 


9xik[x,ek(xH  =  cpxlk(x,0)  +  ek(x)cpxyik(x,0)  + 


cpylk[x,ek(x)]  =  qpylk(x,0)  +  €k(x)cpyyik(x,0)  +  . 


(34) 


Making  use  of  equations  (34),  substituting  equation  (6)  into  equa¬ 
tion  (33),  and  ordering  in  powers  of  b  and  e,  the  following  set  of 
boundary  conditions  results: 


-In  b: 


In  e : 


In  b2: 


cpy10(x,0)  =  -cpy01(x,0)  =  -n(x) 


cpy12(x,0)  =  k'(x)cpx00  =  k'(x)Mo 


cPy11(x,0)  =  -cpy02(x,0)  =  0 


(35a) 


(35b) 


(35c) 


In  eb: 

-cpy13(x,0)  =  k'(x)[9x01(x,0)  +  qpx10(x,0)]  -  k(x)|^yy01(x,0)  +  <pyy10(x,oi] 
=  k'(x)|m(x)  +  cpx10(x,oTj  -  k(x)[p2m'(x)  +  4yy10(x,o7|  (35d) 

In  e2: 

cPylJ+(x,0)  =  k'  (x)cpx12(x,0)  -  k(x)cpyy12(x,0)  (35e) 

Also  in  equations  (35)  the  right-hand  sides  are  known  after  the  lower- 
order  problem  has  been  solved. 
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Characteristic  Condition 

Procedure.-  It  may  he  noted  that  for  the  computation  of  the  basic 
flow  field  both  velocity  components  were  prescribed,  separately,  along 
the  axis,  whereas  in  the  present  problem  only  the  ratio  of  the  two  veloc¬ 
ities  is  fixed  by  the  airfoil  condition.  Thus,  an  additional  condition 
is  required. 

From  the  fundamental  laws  of  supersonic  flow  it  is  evident  that  no 
perturbation  of  the  basic  field  can  exist  upstream  of  the  two  character¬ 
istics  passing  through  the  leading  edge  in  a.  downstream  direction.  Exist 
ence  theorems  for  hyperbolic  differential  equations  require  that  the 
potential  be  continuous  everywhere.  Hence,  the  perturbation  potential 
must  vanish  all  along  the  leading-edge  characteristics. 

To  formulate  this  characteristic  condition  analytically,  the  equa¬ 
tion  of  the  characteristics  has  to  be  known.  Suppose  the  differential 
equation  is  given  by: 


A(x>y;  Vcpy)?>xx  +  2B(x,y;  <px , q>y) q>xy  + 


c(x,y;  cpx,<Py)9yy  +  F(x,y;  Tx,Ty)  =0  (36) 

where  A,  B,  C,  and  F  are  linear  or  nonlinear  functions  of  x,  y, 
and  the  first  derivatives  cpx  and  cpy.  Then  the  equation  of  the  char¬ 
acteristics  is  as  follows: 

A  dy2  -  2B  dx  dy  +  C  dx2  =0  (37) 

so  that 

(38) 


Thus,  the  characteristics  are  seen  to  depend  only  on  the  coefficients 
of  the  second-order  derivatives  in  equation  (36).  In  general,  these 
coefficients  are  functions  not  only  of  x  and  y  but  also  of  the 
desired  function  cp.  Hence,  in  a  rigorous  theory  the  characteristics 
are  not  known  beforehand  but  only  after  the  function  cp  is  determined. 
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By  assuming  a  power-series  solution  for  cp,  this  indeterminancy  is 
removed;  the  nonlinear  equation,  equation  (36),  is  replaced  by  a  set 
of  linear  differential  equations  whose  second-order  derivatives  have 
the  same  constant  coefficients  for  each  equation.  In  this  way,  the 
characteristics  of  each  linear  equation  of  the  set  turn  out  to  be  the 
same  straight  Mach  lines,  x  ±  Py  =  Constant,  determined  by  the  recti¬ 
linear  component  of  the  basic  flow  field,  whereas  the  true  characteris¬ 
tics  are  functions  of  the  entire  flow  field  (including  the  perturbation 
field  cp’  )  and  hence  they  are  curved.  This  is  an  unsatisfactory  state 
of  affairs. 

It  is  possible  to  improve  the  computation  of  the  characteristics 
somewhat  by  considering  once  more  nonlinear  equation  (5)>  before  it  is 
split  up,  and  by  merely  distinguishing  between  the  potential  of  the  basic 
field  cpB,  which  is  known,  and  that  of  the  disturbance  field  cp'  ,  which 
is  looked  for.  If  the  order  of  magnitude  of  the  various  potential  coef¬ 
ficients  (p1^  is,  for  the  moment,  left  out  of  consideration,  curved 
characteristics  can  be  obtained  also  for  linear  differential  equations, 
as  the  following  illustration  may  show. 

Substituting  the  potential  series  (6)  in  the  form: 

<p  =  cp6  +  cp'  (39) 


equation  (5)  becomes: 


Since  cpB  has  been  determined  already  as  a  function  of  x  and  y,  equa¬ 
tion  (40)  may  be  considered  as  a  differential  equation  for  the  -unknown  cp' 
only  and  may  then  be  written  in  the  following  way: 

Acpxx'  +  2Bcpxy'  +  Ccpyy'  +  F(x,y;  cpx’  ,cpy'  )  =  0 
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Here  F  is,  as  before,  a  given  function  of  x,  y,  and  the  first  deriv-. 
atives  cpx'  and  qpy*  .  Each  of  the  coefficients  A,  B,  and  C  consists 

of  one  part  that  depends  only  on  cpB  and  a  second  part  which  is  propor¬ 
tional  to  cpx’  or  cpy'  : 


(42) 


Since  it  is  the  prime  purpose  of  this  analysis  to  determine  the  influence 
of  the  inhomogeneous  field  upon  the  pressure  distribution  around  the  air¬ 
foil  by  means  of  a  linear  theory,  it  is  sufficient  to  approximate  the 
equation  for  the  characteristics,  by  neglecting  in  the  coefficients  A, 

B,  and  C  that  part  which  depends  on  qp*  ,  so  that  now: 


A 

B  [■ 

C 


f^tp6)  «  gx(x,y) 


In  this  way,  the  characteristics  of  the  total  field  around  the  airfoil 
are  approximated  by  the  characteristics  of  the  basic  inhomogeneous  field. . 
On  the  other  hand,  the  lines  represented  by  equation  (38)  which  result 
from  such  a  procedure  will  be  curved.  Thus,  the  equation  for  the  char¬ 
acteristics  becomes: 


£  +  i1ri§=2  -  (<PyB)2]  -  +  2('fX:ETyB)  d-<  <3y  +■ 

£  +  -  Kb)3  -  I*2  =  0  (43) 


where 


cp8  =  Mqx  +  bcp^l  +  b2cp0^ 


(44) 
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Analytic  formulation.-  Substituting  equation  (44)  into  equation  (43) 
and  retaining  only  terms  up  to  first  powers  in  b,  the  coefficients  of 
equation  (43)  become: 


A  = 


7  +  l)Mobcpx01 


+ 


2B  =  -2M0bcpy01 

C  =  [l  -  (7  -  DMobtpx01  .  .  f] 


(45) 


where  cpx®-*-  and  cpy^l  are  functions  of  m(x,y)  and  n(x,y)  given  by 
equations  (24). 

Inserting  equations  (45)  into  equation  (38)  and  keeping  consistently 
only  terms  up  to  first  powers  in  b,  the  equation  of  the  characteristics 
becomes : 


d y  =  MptKPy01  t  \)[l  -  (7  -  l)Mobqpx0^[3g  +  (7  +  PM^0]] 

dX  B2  +  (7  +  l)Mobcpx01 

and  after  expanding  the  right-hand  side  in  powers  of  b: 


-b|lE  + 1  +  -  ^ik01  ±bir  V31  •  •  •  }  (h7> 

The  upper  (or  lower)  sign  is  used  for  the  characteristics  above  (or  below) 
the  airfoil. 

Equation  (47)  may  be  solved  by  an  iteration  process,  in  which  the 
first  approximation  is  obviously 


dy  _  +1 

dx  P 
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This  leads  to  the  Ma.ch  lines: 

y  =  ±x/(3  (48) 

where  the  constant  of  integration  has  been  set  equal  to  zero  so  as  to 
yield  the  leading-edge  characteristics. 

Substituting  equation  (48)  into  the  right-hand  side  of  equation  (47) 
yields  in  the  second  approximation: 


(49) 


From  equations  (24)  substitute  for  the  characteristics  above  (or  below) 
the  airfoil: 


('Oy-tx/p 


=  it&O  ±  ♦  m(0)  X  hMI 

2\_  (3  P  J 

=  |£m(2x)  ♦  +  m(0)  +  2^1] 


>• 


(50) 


Inserting  equations  (50)  into  equation  (49)  and  integrating  yield 
for  the  characteristics: 


y  =  70(x)  +  bt1(x)  .  .  . 


where 


t0(x)  =  ±x/P 


The  upper  (or  lower)  sign  refers  to  the  upper  (or  lower)  surface, 
respectively. 
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Coming  back  to  the  characteristic  condition  indicated  in  the  sec¬ 
tion  "Procedure,"  the  following  requirement  has  to  be  fulfilled:  For 


y  =  i|r0(x)  +  b^1(x)  +  .  .  . 


cp'  =  bcpb®  +  b2cpbb  +  ecpb2  +  ebcpb3  +  e2cpl4  =  0 
More  explicitly: 

bcp10[x;  i|r0(x)  +  btp(x)  .  .  7j  +  b2cp1;L^;  \|rQ(x)  +  bi1(x)  .  .  .1  + 
ecp12[x;  i(r0(x)  +  bi^x)  .  .  7|  +  ebcp13[lc;  ^(x)  +  b^x)  .  .  7j  + 

40(x)  +  b4]_(x)  •  •  n  =  ° 

Expanding  qp'  about  y  =  4q(x)  in  Powers  of  b  yields: 

bcp10(x^0)  +  b2|t1(x)cpy10(x^0)  +  ^(x,^)]  +  ecp12(x,t0)  + 

€b  ^1(x)cpy12(x^0)  +  cp13(x,i0)]  +  e2cpl4(x,^0)  =  0 

Splitting  equation  (53)  in  powers  of  b  and  e  furnishes  for  each 
potential  coefficient  cp1^  one  characteristic  condition: 


(52) 


(53) 


In  b: 


In  e : 


In  b2: 


cP10(x^o)  =  0 


<P12(x^o)  =  0 


(54a) 


(54b) 


<P11(X,t0)  =  -tl(X)(Py10(X^4'o) 


(54c 
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In  eb: 


In 


q)13(x^0)  =  -t1(x)cpy12(x^0) 


(54a) 


(54e) 


Whenever  the  relation  is  inhomogeneous,  the  unknown  potential  is  expressed 
in  terms  of  known  lower-order  potentials. 

Each  of  the  boundary- value  problems  consisting  of  a  differential 
equation,  a  boundary  condition,  and  a  characteristic  condition  has  now 
been  formulated  and  will  be  solved  in  the  appropriate  order. 


COMPUTATION  OF  PERTURBATION  POTENT IAIS 


In  the  following  discussion  (see  fig.  l)  k(x)  denotes  the  ordinate 
of  the  profile  above  (or  below)  the  x-axis,  which  is  the  direction  of  the 
relative  wind  Mq.  To  distinguish  between  the  upper  and  lower  surfaces, 

ky  and  kp  are  used;  k(l)  is  zero  or  unequal  to  zero  depending  on 

whether  the  airfoil  is  at  zero  or  at  finite  angle  of  attack.  If  there 
are  two  signs  in  front  of  a  term,  the  upper  sign  refers  to  the  upper, 
and  the  lower  sign,  to  the  lower  surface. 


Computation  of  cp-*-® 

That  part  of  the  disturbance  potential  which  is  proportional  to  b, 
namely,  qP-O,  obeys  the  following  equations:  The  differential  equation: 


P^xx10  -  <Pyy10  =  0 

the  boundary  condition: 

9y10(x,0)  =  -n(x) 

and  the  characteristic  condition: 

cp10(x,±x/p)  =  0 


(55) 
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To  obtain  a  solution,  it  is  convenient  to  introduce  the  characteristic 
coordinates : 


I  =  x  +  3y 
T]  =  x  -  (3y 


Then 


and 


=  <PS  +  ^ 

9y  =  0(q>|  -  <Pr,) 


(56) 


B^cp  -  m 
H  Yxx  uyy 


(57) 


As  the  characteristic  condition  is  different  on  the  upper  and  lower  sur¬ 
faces,  the  two  surfaces  will  be  treated  separately. 

Upper  surface.-  In  characteristic  coordinates  equations  (55)  maybe 
written  as  follows  for  the  upper  surface:  The  differential  equation 
becomes 


4pV  10  =  0 

the  boundary  condition, 

"(V0  -  ■  -<*> 


(58) 


and  the  characteristic  condition, 

cp10(|,0)  =  0 

Integration  of  the  differential  equation  yields: 


<p10(M)  =  F(|)  +  G(t]  ) 


(59) 


where  the  arbitrary  functions  F  and  G  are  determined  by  the  boundary 
condition  and  the  characteristic  condition. 
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From  the  characteristic  condition  there  follows: 

=  F' (!)  =0  for  any  | 


(60) 


Then  the  boundary  condition  requires  that 


pfF'm  -  G'(ti[|^=x  =  p[f’  (x)  -  G'(xj]  =  -n(x)  (6l) 


Taking  into  account  equation  (60),  it  follows  that: 

G’ (x)  =  i  n(x) 

r 

G'  (rj)  =  -  n(q) 


(62) 


Therefore: 


9x10U,y)  =  |  n(x  -  py) 
Ty10(x,y)  =  -n(x  -  py) 


(63) 


To  obtain  the  pressures  on  the  airfoil  surface  y  =  ek(x),  equations  (24) 
may  be  expanded:  '  ' 


=  cpx10(x,0)  +  eku(x)cpxy10(x,0) 

=  ~  n(x)  -  ekylxjn1 (x)  +  .  .  . 


>  m 


Ty10 ^6ku(x)]  =  -n(x)  +  e3ku(x)n>  (x)  +  .  .  . 


same: 


Lower  surface.-  On  the  lower  surface,  the  boundary  condition  is  the 


Kn10  -  tft]10) 


-  -  -w 
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but  the  characteristic  condition  differs: 


From  equation  (65)  there  follows:  cp^(0,q)  =  G' (t|)  =  0,  for  any  T) . 
Then,  from  the  boundary  condition  there  follows: 


f3cp^  =  PF'(|)|=X  =  -n(x) 


so  that 


*6(S) 

Using  equations  (4-2) 


(66) 


and 


q>x10(x^y)  =  n(x  +  Py) 


m  10/  \ 

9y  U>y) 


-n(x  +  Py) 


(67) 


On  the  airfoil  surface: 


ekL(x)n' (x)  .  .  . 


Cpy10  ^ekL(x[]  =  -n(x)  -  epkL(x)n’  (x)  .  .  . 


(68) 
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Computation  of  cpH 

For  the  potential  cpH  which  is  proportional  to  b2,  the  boundary- 
value  problem  ca.n  be  formulated  by  means  of  equations  (32c),  (35c),  ' 
and  (54c):  The  differential  equation  is 

-  V11  ■  -”o  {lir  +  i)  +  (y  -  l)P^](9xx01+  cPxx10)(cRx0;L  +  9x10)  + 

2(<Py01  +  9y10)(<Pxy01  +  V10)}  +  V02  -  P2<PXX02 

the  boundary  condition, 

cPy11(x,0)  =  0 

and  the  characteristic  condition, 
cp1^(x,±x/(3)  =  -^1(x)cpy10(x^±x/p) 


Upper  surface.-  Inserting  the  values  for  cp01,  cp1^,  and 

op  p  op 

<Pyy  -  P  9xx  from  equations  (24),  (64),  and  (30),  for  the  upper 

surface  problem  (69),  in  characteristic  coordinates,  is  as  follows: 
The  differential  equation  becomes 


4p2cp. 


11  Mo 


In 


Hr  +  I)(e2  +  i)]^U)2hiU  m. (Tjjabil]  x 
[Tr  +  l)  +  (r  -  3)pf|  +  hUJ  + 


n(rj) 

P 


(70) 


the  boundary  condition. 


h11  -  V1) 


|=T]=X 


=  0 


(7D 


and  the  characteristic  condition. 


<P1:LU,0)  =  i|f1(|,0)n(0) 


(72) 
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The  rest  of  the  procedure  may  be  indicated  here  briefly,  the  details 
being  carried  out  in  appendix  A. 

Upon  integrating  equation  (70)  with  regard  to  |  and  tj,  respectively: 
cd^11  =  G1(m,n)  +  P-^U) 

<PT111  =  02(m,n)  +  P2'  (q) 

where  G^  and  G2  are  known  functions  of  m(  t ) ,  m(q),  n(|),  and 
n('q )  and  P-^  and  P2  are  unknown  integration  constants. 

Now,  characteristic  condition  (72),  upon  differentiating  with  regard 
to 


(73) 


9|1;L(^0)  =  t]^U,0)n(0)  (74) 


determines  Pj'^)  by  means  of  equations  (73). 

Then  boundary  condition  (71),  after  insertion  of  equations  (73), 
fixes  P2'  (q). 

With  P^  and  P2  determined,  the  two  velocity  components  cpxll(x,0) 
and  cpy-^(x,0)  that  are  necessary  for  the  computation  of  the  pressure  on 
the  airfoil  surface  can  be  determined  by  means  of  equations  (73)  and  (56). 

Lower  surface.-  For  the  lower  surface  an  analogous  procedure  can  be 
carried  through;  note  that  the  characteristic  condition  is  now: 


cPt^11(°,t1  )  =  4f1(0,q  )n(0) 
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The  resulting  velocity  components  on  either  airfoil  surface  are 
as  follows: 


-  3)Mo2  +  3'*^ +  E 


n'(x) 


~>x 


gTV  ) 


dV 


+  m’  (x)  + 


[  +  2^y  -  1)M02  +  g*l*>l*I  ■+ 


\(7  +  DM0^)- 


(0)‘ 


9v11(xJ0)  =  0 


y  (75) 


Computation  of  cp12 

That  part  of  the  disturbance  potential  which  is  proportional  to  6, 

12 

cp  ,  must  be  expected  to  equal  the  expression  obtained  in  the  Acker et ' 
theory. 

Hence,  the  velocity  components  above  (or  below)  the  airfoil  are: 


?x12  =  +y  ku(L),(x  +  Py) 


V2  =  Mo%(L)’(x  +  py) 


(76) 


By  expanding  about  y  =  0,  the  velocities  on  the  upper  (or  lower) 
airfoil  surface  are: 


M, 


q>x12&^(x3  =?jVl)'W  +  MoekU(L)(x)ku(L)”(x) 


cpy12g,sk(x3  =M0ku(L)*(x)  ?  M0P6ku(L)(x)ku(L)*'(x) 


(77) 
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Computation  of  cp^3 

The  potential  cp^-3,  which  is  proportional  to  eb,  may  be  considered 
to  originate  from  some  kind  of  interaction  between  the  inhomogeneity  of 
the  basic  field  a,nd  the  airfoil  disturbance.  The  boundary-value  problem 
for  cp!3  can  be  formulated  b;y  means  of  equations  (32d),  (35b),  and  (54d) 
The  differential  equation  is 


32Txx13  -  Tyy13  =  -M.  J 


^[7  +  1  +  (7  -  1)p5  9x12(cpxx01  +  9XX10)  + 


^Pxx12^01  +  <P:  10 


)]  +  2  l?y12(< 

i 


'Pxy01  +  ^xy10)  + 


CPxy12(<Py01  +  CPy10 
the  boundary  condition, 

Cpy13(x,0)  =  (x)[m(x)  +  ,x10(x,03  " 

kU(L)^x)  [P2m' +  cPyy1°(x^c0] 
and  the  characteristic  condition, 

cp13(x,±x/p)  =  -i1(x)cpy12(x,±x/p) 


W78) 


Upper  surface.-  The  treatment  of  the  upper  surface  is  as  follows: 

Inserting  for  cp-^  and  cp^2  expressions  (63)  and  (77),  problem  (78) 
written  in  characteristic  coordinates  becomes: 

The  differential  equation: 


^3  - 

in 


Mo 

2P 


E?  +  1)  +  (7  -  3)p5  jy  (n)[m’(l)  +  — +  ky"  (n)jm(l)  +  iUl  j. 
E y  +  D(P2  +  !)]  •jjcu' (ri)jm'  (n)  +  +  WCn^n)  + 


+ 


(79) 


NACA  TN  2849 


31 


* 


The  boundary  condition: 


p(^13 '  b13)t=1 


V U) 


m*  (x) 


n’(x) 

P  _ 


(80) 


The  characteristic  condition: 

.  cp13(M)  =  -^1(|,0)Moku’(0)  (=  0) 


(81) 


The  rest  of  the  procedure  may  be  indicated  here  briefly,  the  details 
being  carried  out  in  appendix  A. 

Upon  integrating  equation  (79)  with  regard  to  |  and  t),  respectively: 

cp|13  =  F1(m,n,k)  +  J-l’ (I)  (82a) 

and 

cpq13  =  F2(m,n,k)  +  Jg' (q)  (82b) 

where  Fp  and  F2  are  known  functions  of  k(|),  mU),  n(i),  k(q), 

m(q),  and  n(q)  and  J-^  and  Jg  are  unknown  integration  constants. 

Now  characteristic  condition  (8l) 


9^3(1, 0)  =  -ti|a,0)Moku*(0) 


determines  J^' (| ) 

vanishes  for  k(0) 
Also  equation  (8l) 


by  means  of  equations  (82).  The  function  Fj_ 

=  0  and  k1  (0)  =  0,  a,s  can  be  seen  from  appendix  A. 
vanishes  in  that  case.  Then  it  follows  simply  that: 


Jp*  U)  =  0 

This  simplification  suggests  making  the  assumption  k' (0)  =  0,  which  is 
an  artifice  frequently  used  in  linearized  theories.  The  error  introduced 


4. 
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by  this  deviation  from  the  true  contour  is  known  to  yield  small  inac¬ 
curacies  in  the  flow  field  in  a  very  thin  strip  adjacent  to  the  two 
characteristics  passing  through  the  leading  edge.  Therefore  the  sim¬ 
plification  is  permissible. 

Boundary  condition  (80),  after  insertion  of  equations  (82),  fixes 
Jo'  (q).  With  Jj  and  J2  fixed,  the  two  velocity  components  cpx13(x,0) 

a.nd  cpy-*-3(x,0)  that  are  necessary  for  the  computation  of  the  pressure  on 
the  airfoil  surface  can  be  determined  by  means  of  equations  (82)  and  (56). 


Lower  surface.-  For  the  lower  surface  an  a.nalogous  procedure  can 
be  carried  through.  Note  tha.t  now  the  characteristic  condition  is 
q>r|(0,T))  =  0  and  that  different  expressions  for  k(x),  cplO,  and  cp-^ 
have  to  be  inserted. 


The  resulting  velocity  components  on  the  airfoil  surface  are: 


13(x,0)  =  ±^£  JlkM  2(7  -  1)  +  8  - 


8p2 


*oc. 


ku(L)’(x)[®(x) 4  ^r]  ■ 

-3)  +  8  +  2fi|ku(L)U)[5(x)  ±  j.  (83) 

cpy13(x,0)  =  ku(L)'(x)j^i(x)  ±  ^Llj  _  32ku(L)(x)  m'  (x)  ±  n  ^  J 

Computation  of 

p 

That  part  of  the  disturbance  potential  which  is  proportional  to  e  , 
that  is,  cpl\  obeys  equations  (32e),  (35e),  and  (5^e):  The  differential 
equation  is 


W4-  14  - 


Tyy1^  =  "Mo[(7 +l)cpxx12cPx12+  1)CP: 


,  12  +  2qj  12_  12 

X  tyy  T  TXy 


the  boundary  condition, 


cpyll<-(x,0)  =  Tx12^0)^!)'  (x)  -  9yy12(x,0)ku(L)(x) 
and  the  characteristic  condition, 

<P1^(x,±x/p)  =  0 


ym 
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Upper  surface.-  For  the  upper  surface,  after  inserting  for  cp12 
expressions  (77)>  problem  (84),  written  in  characteristic  coordinates, 
is  as  follows:  The  differential  equation  becomes 

4p2cp|7)l4  =  -  p-|^7  +  l)V(3)V’(q)]  (85) 

the  boundary  condition, 

-  (Pnl4)l=T1=x  =  -  ijr  V2W  +  PM0ku(x)ku««(x)  (86) 
and  the  characteristic  condition, 

<pl4U,0)  =  0  (87) 

Integrating  equation  (85)  with  respect  to  |  and  q,  respectively: 


“  "  +  Mi'^' 


14  ^5 


^-£7  +  l)V(q)ku-(q)r|  +  M2-(q) 


where  M-^  and  Mg  have  to  be  determined  from  the  boundary  condition 
and  the  characteristic  condition. 

From  characteristic  condition  (87)  there  follows: 

qptl4U,o)  =  0  (89) 

which,  by  means  of  equations  (88),  fixes  M-j/U)  =  0  if  again  the  assump¬ 
tion  k'  (0)  =0  is  made,  as  in  the  computation  of  cp!3# 

Then  boundary  condition  (86),  after  insertion  of  equations  (88) 
fixes  M2‘ (q). 

With  and  M2  fixed  the  two  velocity  components  on  the  airfoil 

surface  can  be  determined  by  means  of  equations  (88)  and  (56). 
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Lower  surfa.ce.-  After  analogous  considerations  for  the  lower  surfa.ce 
ha.ve  been  carried  out,  the  velocity  components  on  the  two  surfaces  are: 

cpxl^(x,0)  =  -i.  jj(7  +  1)M0J+  -  4(^2  -  l)\  ku(L),2(x)  + 
^pi4lcU(L)^x)kU(L)" 

cPy1^(x^°)  =  +"j^  kU(L)'2^x^  ±  PH3kU(L)(x^kU(L)"  ^x^ 

FORCES  AND  MOMENTS  ON  AIRFOIL 
Complete  Velocity  Field  oh  Airfoil  Surface 

To  obtain  the  pressure  distribution,  the  squares  of  the  velocity 
components  of  the  whole  field  are  needed.  By  adding  all  of  the  com¬ 
ponents  (up  to  the  second  order  in  b  and  e),  there  result: 


Substituting  the  expressions  derived  in  the  preceding  sections: 


2-2  2  2 
cpx  =  MQ  +  bA1  +  eAg  +  b  A^  +  ebA^  +  e 


cpy2  =  e2M02ku(L)l2(x) 


(92) 


Pressure  Distribution  on  Airfoil  Surface 

The  pressure  distribution  can  be  obtained  from  the  velocity  distri¬ 
bution  by  means  of  the  energy  equation  (reference  2)  which,,  in  dimension¬ 
less  form,  is  as  follows: 
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Here  p  ,  pQ,  and  Mq  are  taken  for  the  homogeneous  flow  field  and  p, 

p,  q,  and  c  are  local  quantities.  By  use  of  the  isentropic  p  -  p 
relation,  equation  (94)  becomes: 


and  after  substituting  p  =  pQ 

relative  pressure  increment  — 

Po 


+  Ap  and  expanding  in  powers  of  the 
■,  which  is  small  compared  with  1: 


Upon  substituting  equations  (92)  for  cpx2  and  cpy2,  equation  (96) 
may  be  solved  in  successive  approximations.  In  the  first  approximation 
and  second-order  quantities  in  the  velocities  are  neglected. 

VVo) 

Then: 


Ap 


Poco 


-  =  -bMn 

2  0 


m 


(x)  ± 


n(x) 


V 


kU(L)’(x) 


(97) 


Upon  substituting  from  equation  (97)  for 


and  keeping 


second-order  quantities  in  the  velocities  throughout  equation  (96),  the 
second  approximation  yields: 


€P1  +  bP2  +  b2P3 


+  ebP^  +  6 
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Here 


As  shown  in  appendix  B,  equations  (98)  for  the  pressure  distribution 
on  the  airfoil  surface  can  be  checked  in  a  limiting  case;  in  which  com¬ 
plete  agreement  of  equations  (98)  with  a,n  expression  derived  by  Ferri 
(reference  2)  is  obtained. 

Drag  Force 

In  order  to  obtain  the  drag  force,  the  component  of  the  pressure 
which  is  parallel  to  the  flight  direction  x  has  to  be  integrated  over 


38 
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the  upper  and  lower  surfaces.  Denoting,  as  in  figure  1,  by  e©u(p)  "the 

angle  subtended  by  the  local  tangent 
direction,  the  drag  becomes: 

V  =  J sin  e  6 y  dSy 

where  ds^  and  ds^  are  the  line  elements  of  the  airfoil  surface. 

Keeping  only  terms  up  to  third  orders  in  e,  equation  (99)  becomes: 


of  the  surface  with  the  flight 


Ap  sin  e 


6L  dSL 


(99) 


D 


f  ApyyU) 


dx 


Ap^ek^' (x)  dx 


Hence,  the  drag  coefficient  is: 


%Qg 

|poMo2! 


Ap 


2^oco 


2 


(100) 


For  (Ap)TT  and  (Ap)T  expressions  (98)  have  to  be  substituted  in 
U  -L 

equation  (100). 


In  the  case  of  a  symmetrical  airfoil  at  zero  angle  of  attack 
ky(x)  =  -kL(x)  =  k(x).  Then: 


(cd) 


=  e  D 


6=0 


+  1o€D2  "** 
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Here: 


P  'Jq 


B2=7T  m’ (x)k(x)  ^ 
Mo  J0 


D3  =  — f-o 
Mn2p2 


Or  -  1)0+  2]  f  ^  00  dx  * 

Jo  P2 


[(7  -  3)M04  +  i+.MQ2| jT1  HIM  k*  (x) 


°-%  P2 


2^M02  -  l)2  /  m2(x)k'  (x)  dx 

J0 


n(x* )  dx' 


dx  + 


D4  = - =■ 

MqPj 


\(y +  DMo4  -  k{%2  - !)] 


2  /  m(x)k,2(x)  dx  + 

J0 


m'  (x)k(x)k'  (x)  dx 


D5  =^1 


(7  +  l)M01+  -  4(mo2  -  lj  J  k’ 3(x)  dx 


(101) 


In  the  case  of  a  symmetrical  airfoil  at  a  finite  angle  of  attack  e0 
the  decomposition  of  ky' (x)  and  ky'(x),  suggested  in  the  section  on 
airfoil  geometry  (appendix  C),  has  to  be  carried  out.  In  other  words: 
Replace  ky’  in  (Ap)y  and  ky'  in  (Ap)y  by: 


and 


ky’ (x)  =  k(x)  -  e 


ky* (x)  =  -k(x)  -  6 


(102) 
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where  now  k(x)  sj  k(x)  is  the  customary  profile  function,  for  which 
k'(l)  =  0. 

Then,  upon  substituting  the  modified  expressions  (98)  into  equa¬ 
tion  (100),  the  following  generalizations  of  the  coefficients  of  Cp 
have  to  he  noted: 


6C 
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1+1 

The  physical  meaning  of  the  various  terms  in  equations  (101) 
or  (103)  may  be  readily  interpreted:  D1  (or  D1' )  is  the  wave  drag 

of  the  airfoil  well-known  from  the  Ackeret  theory;  (or  De>' )  is 

the  second-order  drag  known  from  the  Busemann  approximation  (refer¬ 
ence  3);  Dg  represents  the  horizontal  buoyancy  acting  on  the  airfoil; 

may  be  denoted  "second-order  buoyancy";  and  is  the  term  charac¬ 

terizing  the  interaction  between  the  airfoil  field  and  the  given  pres¬ 
sure  gradient  of  the  wind  tunnel. 

Insofar  as  beDg  is  concerned,  it  is  to  be  noted  that  the  integral 
is  taken  over  the  product  of  the  local  volume  element  ek(x)  dx  and  the 

ii  .  . 

local  pres sure -gradient  coefficient  b —  m' (x).  Since  the  pressure  gra- 

M0 

dient  is  varying  from  point  to  point,  this  expression  represents,  indeed, 
the  generalization  of  the  classical  result  which  is  thus  seen  to  hold  not 
only  for  incompressible  but  also  for  supersonic  flow. 

However,  the  other  correction  term  for  the  drag,  which  can  be  derived 
for  incompressible  flow  (see  "Introduction"),  does  no  longer  appear  in 
the  supersonic  range:  The  simple  concept  of  the  "apparent  mass"  is  to  be 
replaced  by  more  complicated  expressions. 

Lift  Force 

In  order  to  obtain  the  lift,  the  component  of  the  pressure  which  is 
perpendicular  to  the  flight  direction  has  to  be  integrated  over  the  upper 
and  lower  surfaces  of  the  airfoil.  Thus: 


L  =  -J  (Ap)jj  cos  e  dSy  + 
which  up  to  and  including  second  orders  in  g  becomes: 


(Ap)  L  cos  e  6 k  ds^ 


[(AP)U  -  (Ap)J  dx 


L  x  -Z 


(104) 
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The  lift  coefficient  is: 


CL  = 


Mn2j0 


Ap 


Ap 


&oco2)jj  \hcoVL 


dx 


(105) 


For  (Ap)„  and  (Ap)^  expressions  (98)  have  to  be  substituted  in 
equation  (105). 

In  the  case  of  a,  symmetrical  airfoil  at  zero  angle  of  attack 
ky(x)  =  -k  (x)  =  k(x).  Then: 


o 

bl^  +  b  L2  +  ebL^ 


Here: 


L 


1 


-  3)Mr 


4m, 


'0 


n(x’  ) 
P 


dx  + 


2 


+ 


DM,5  -  2(m02  - 


m(x) 


*7 


M106) 


Lo  = 


Mnp3 


E 


(7  +  DMr 


n(x) 

P 


k’ (x)  dx  + 


Note  that  equations  (106)  vanish  for  m 


J 

n  =  0,  as  they  should. 
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In  the  case  of  a  symmetrical  airfoil  at  a  finite  angle  of  attack  6 , 
after  making  the  same  modifications  as  in  the  case  of  the  drag,  the  fol¬ 
lowing  changes  in  the  coefficients  Lj_  have  to  be  made: 


L 


1 


L 


1 


*  —  Ij^ 


L3  " 


MqP: 


[(7  +  1)M04  -  4(^2  -  lj] 


2  f  m(x)  dx  + 
J0 


(107) 


dx 


Furthermore,  a  coefficient 


appears . 


e 


4e 

P 


The  coefficients  L^'  and  Lr>'  represent  the  lift  expressions 
known  from  the  work  of  Ackeret  and  Busemann  (see  reference  3), 

and  L2  a.re  vertical  forces  due  mainly  to  the  vertical  pressure  gradient, 
and  is  the  interaction  term. 


Pitching  Moment 


The  moment  about  the  leading  edge  is  obtained  by  integrating  over 
the  moments  of  the  differential  lift  forces  dL  acting  at  a  distance  x. 
Therefore  the  moment  coefficient  becomes: 


M/c02 

—  P  M  ^2 
200  L 


Ap 


Mq2 


J0  Ipc  2 1 

\2hb  o  ) 


\x  dx 


The  integral  is  to  be  taken  over  the  upper  as  well  as  the  lower  surface. 


Since  an  upward  lift  which  is  to  be  considered  positive  produces  a 
negative  (diving)  moment,  the  contribution  of  the  upper  surface  is: 
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The  contribution  of  the  lower  surface  is: 

U 


(cm)l  =  - 


1 


Ap 


x  dx 


MO  0  \oPnC 


2^o  o  /l 


Hence : 


CM  "  w  2 


1_ 

V 


PI 


Jo 


Ap 


Ap 


■  TTPnC  2 


Ip  C  2 


2H°  o  /u  \2  o  o  /L 


x  dx 


In  the  case  of  a  symmetrical  airfoil  at  zero  angle  of  attack 
kyCx)  =  -kL(x)  =  k(x) .  Then: 


bM1  +  b2M2 


+  €  IdM^ 


Here: 


(108) 


>(109) 
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In  the  case  of  a  symmetrical  airfoil  at  a  finite  angle  of  attack  8 
after  making  the  same  replacements  as  before: 


Furthermore,  two  additional  coefficients  appear: 


and 


r _ 

— . 

e2M5-  =  e2  2S  < 

5  I 

j|p  +  1)M04  -  ■ 

4mo2  - 1 

)  /  xk'  (x)  dx 

_jj0 

The  coefficients  M^'  and  M^’  represent  the  moments  known  from 
the  Ackeret  and  Busemann  calculations,  and  M2  are  the  moment  cor¬ 
rections  due  to  the  pressure  gradients  of  the  wind  tunnel,  and  is 

the  interaction  term. 
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Application  of  Theory 

As  a  simple  example  of  the  foregoing  general  considerations,  the 
force  corrections  will  be  given  for  a  symmetrical  airfoil  placed  at 
zero  angle  of  attack  into  a  tunnel  field,  for  which 


m(x)  =  n(x) 


x 


^i.e.,  for  which  the  pressure  gradient  bpUm' (x)-y-  is  linear  in  xj: 


=  e2  —  /  [k1  (x)J^  dx  +  eb  —  /  k(x)  dx  + 


3  J0 


Mo  J0 


e3  ~[(7  +  1)M04  -  4(mo2  -  1  )\J  (k'U)]3  dx  - 
eb2  ^-[m02(37  -  17)  +  i*(mo1*  *  **3  i  xk(x)  dx  - 


e2b  -JLJfty  +  l)^4  -  4(mq2  -  ljl  f  xk'2(x)  dx 
Mop-3  L  -1  J0 


>  (HD 


2  t,2  1 


Ct  =  b  —  —  b 

L  MoP  Buffi 


M, 


4(57  +  l)  -  4(^2  -  2)]  - 


eb  -—-[(7  +  1)Mq4  -  4(mq2  -  1 )\J  k(x)  dx 


MqP 


Cm  -  -b 


h +  b2  ♦ 1)  -  “ta2  -  23 


3M03  Sm^3 


eb  -1_|mo4(7  +  1)  -  4(mo2  -  lj]  f  xk(x)  dx 


MqP 
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For  a  doubly  symmetric  diamond  profile  of  vertex  angle  e  these 
formulas  reduce  to: 


Numerical  results  obtained  by  means  of  equations  (112)  are  plotted 
in  figures  2  to  9  for  representative  values  of  b  and  e.  In  figure  2 

the  components  of  C-q  which  are  proportional  to  e2,  eb,  e2b,  and  eb2 

are  graphed  separately,  and  in  figure  3  the  total  CD  is  shown  as  a 
function  of  the  Mach  number  M0.  Similarly,  the  components  of  and 

CM  which  are  proportional  to  b,  eb, .and  b2  are  plotted  separately 
in  figures  4  and  6,  respectively,  while  the  total  C-^  and  CM  are  shown 

as  functions  of  Mq  in  figures  5  and  7?  respectively.  For  the  case  of 

a  finite,  small  angle  of  attack  eF  the  behavior  of  the  lift  coefficient 
and  its  components  is  demonstrated  in  figures  8  and  9. 
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Th®  relative  importance  of  the  various  terms  may  be  seen  from  the 
Plots;  .for  example,  for  the  drag  the  first-order  horizontal  buoyancy  is 
the  major  correction  to  be  applied  and  is  on  the  order  of  25  percent  of 
the  true  drag.  Next  in  importance  is  the  interaction  term  which  amounts 

to°  2SS  1  percent  *  ^  ^  The  ^ncy  amounts 

The„Jame  r®lative  magnitude  of  the  correction  terms  exists  in  the 
figure  8fflCient'  lf  the  0536  °f  an  angle  °f  attack  is  considered  as  in 

For  the  lift  and  moment  at  zero  angle  of  attack  the  term  proportional 
small  compared”^ th^it?  ^  lErgeSt  C°rreCti0n'  “hlle  remainder  are 
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APPENDIX  A 

COMPUTATION  OF  PERTURBATION  POTENTIALS 
qp11,  qp!3,  AND  cp1^ 

Computation  of  <pll 

The  details  of  the  computation  of  the  potential  mil 
follows: 

Upper  surface.-  For  the  upper  surface  the  derivatives 
tions  (73)  are  given  hy: 


=  +  1)(1  +  + 

8p^  \  p 


[(7  +  1)  +  (r  -  3)P0-  m'(l)  + 


E 


m(6)  + 


n(!) 

]n(il) 

P  _ 

J  3 

r  +  Pl'(S) 


11  _  Mo 
8p: 


Mr*  Dd*  m'(n)2hl] 


U  +  1  +  (7  -  3)p£] "  m(|)  + 


n*  (tQ 

3 


+ 


are  given  as 
in  equa- 


dq’  + 


(Ala) 


I  + 


(Alb) 
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Equation  (74)  in  explicit  form  is  as  follows: 


<P,U<6,0)  = 

1  8p3 


( 7  +  l)(l  +  P2) 


m 


(0) 


n(0) 

P  J 


+  l  +  (7  -  3)P^] 


m1 


(0  + 


n(  I ) 


-n(O) 


(A2) 


The  constants  of  integration  are  found  to  be: 


Pp’U)  =  +  Dd  +  p2)^^  + 

8p‘ 


[7  +  1  +  (7  -  3)P2] 


m' 


(£)  + 


n’(l) 


n(r|’  ) 


dq ' 


(A3) 


p2’ (h)  =  ~p!  ( 7  +  i)d  +  p^)<! 

8p^ 


n2(0)  m(r|)n(r|) 


m(T|  )-n-,i-3A  +  m' 


>  + 


(a4) 
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Lover  surface.-  For  the  lower  surface  the  derivatives  corresponding 
to  those  in  equations  (73)  are: 


^l11  =  ^2'^7  +  1^1  + 


7  + 


7  + 1  +  (7  -  3)Pg 


m(n) 


n(ri ) 

n(6) 

P  _j 

P 

n'(0  r n 


m(V)  -  ^hlil  dV 


^  +  Pi'd) 


(A5a) 


vx  =  ^((7  +  1)(1  +  + 


11  _ 


[7  +  1  +  (7  -  3)P£ 


m(ri) 


n(q 

find) 

p 

P 

m 


(n) 


-n2¥\P*V-*' 


(A5b) 


The  characteristic  condition  in  explicit  form  is  as  follows: 

Mo 


9T)11(0^)  =  ^3  <j(7  +  1)(1  +  32) 


m(0)  + 


n(0) 

P  _ 


[7  +  1  +  (7  -  3)pgjm(ri)  -  > 


n(  0)  (A  6) 
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In  this  case  the  integration  constants  come  out  to  be: 

Fj/U)  =  -  %f(7  +  1)(1  +  32)Ijm(i)H!lii  +  m’  (| )~~j  I  - 
8p2\  U-  p  13  -I 

£1°)  +  m(5)n(il]j  + 

[7  +  1  +  (7  -  3)3^]  p" -  J  ^  [m(ri •  )  -  -n.^.  ^  dn * 

t",  -  mi:  *?■ »■}) 

p2'(i)  =  -  *  1) <  1  -  P2>[-  ^  * 

\y  +  1  +  (r  -  3)32  jm'  (n)  -  ”  pT| °^-i  al'T 

Computation  of  cp^-3 

The  procedure  for  computing  the  potential  cp!3  is  as  follows.  For 
the  upper  surface  the  derivatives  in  equations  (82)  are  given  by: 


Putting  k(0)  =  k' (0)  =0,  equation  (A9a),  together  with  equa¬ 
tion  (8l),  yields  J^' ( I)  =  0.  On  the  other  hand,  substituting  equa¬ 
tions  (A9a)  and  (A9b)  into  equation  (80)  yields: 


The  lower  surface  can  be  treated  analogously. 
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Computation  of 


In  computing  the  potential  cpl4  the  upper  surface  is  treated  as 
follows.  Putting  in  the  first  of  equations  (88)  k(o)  =  k*  (0)  =  0,  this 
equation,  together  with  differentiated  equation  (87),  yields  M-^'  (£  )  -  0. 

On  the  other  hand,  substituting  equations  (88)  into  equation  (86)  yields: 


Mo’  (tl)  = 


Mo' 


^  +  1  -  2(n)  -  (7  +  i)^'  (n)^"  (ti) 


~^r  ku(q)kuM  (q) 

V 


(All) 


The  lower  surface  can  be  treated  analogously. 
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APPENDIX  B 


CHECK  OF  EQUATIONS  (98)  FOR  PRESSURE 
DISTRIBUTION  IN  A  LIMITING  CASE 


Ferri  (reference  2)  represents  the  pressure  distribution  over  an 
arbitrary  surface,  which  is  exposed  to  a  uniform  stream  of  Mach  num¬ 
ber  M-|_,  as  a  power  series  in  terms  of  the  local  angle  of  attack  q  . 
Thus , 


where  a.^  and  8-2  are  given  functions  of  M-^,  p  is  the  local  pressure^ 
and  pp  the  free-stream  pressure.1 

One  may  note  that  Ferri' s  result  corresponds  to  a  limiting  case  of 
the  general  nonuniform  field,  namely,  to  the  case  in  which  the  basic 
flow  is  characterized  along  the  x-axis  by  an  x-component  MQ  +  bm  and 

a  y-component  bn,  where  both  m  and  n  are  constant. 

The  agreement  of  the  present  result  with  Ferri' s  can  be  seen  as 
follows.  Expressions  (24)  for  cpx01  and  cpy01  indicate  that  if  m(x) 
and  n(x)  are  constant  along  the  axis,  they  are  also  constant  in  the 
whole  xy-plane.  The  same . is  no  longer  true  in  order  b^,  because  the 

differential  equation  for  qp >02  is  inhomogeneous.  However,  up  to 
order  b,  the  basic  field  may  be  considered  uniform  in  this  limiting 


The  Mach  number  of  this  uniform  field  is: 


M1  =  \j(M o  +  bm)2  +  (bn)2 

The  angle  of  the  flow  direction  with  the  local  surface  is 


e0'  + 


Mq  +  bm 


LNote  that  Ferri' s  coefficient  is  here  multiplied  with  M±2. 
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where  e@'  and  r- — - — —  are  the  angles  subtended  by  the  x-axis  with 
M0  +  Dm 

the  local  surface  and  with  the  stream  direction,  respectively. 


In  applying  equation  (Bl)  it  should  be  noted  that  Ferri's  pressure 

coefficient  is  taken  with  respect  to  M-^,  while  the  coefficient  - 2. 

p^c  2 
Ko  o 

in  this  report  is  based  upon  Mq,  which  is  related  to  by  equa¬ 

tion  (B2).  The  connection  between  the  two  pressure  coefficients  is 
given  by  the  relation 


P  -  Pp 

p 

Plcl 


+ 


Pq  ~  Pi 
? 

plcl 


(B4) 


where  — - — ,  representing  the  pressure  difference  between  points 

plcl2 

where  the  Mach  numbers  are  Mq  and  M-j_,  as  well  as  P0/Pi  and  cQ^c1 
may  be  obtained  by  the  isentropic  energy  relation. 

Substituting  Ferri's  result  for  p  -  p-^,  with  and  q  taken 

from  equations  (B2)  and  (B3)  and  the  expression  for  pQ  -  p^  obtained 

from  the  energy  relation,  equation  (B4)  may  be  expanded  in  powers 
of  b  and  yields: 


'2  ~ yj8-  +  -  ‘‘(“o2  -  - 


€b  2(3^^  +  1)M°4  "  "  1ll  (m  *  f)  <B5) 

2 

which  includes  terms  of  order  e,  b,  eb ,  and  e  but  not  terms  of 
order  b2. 
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On  the  other  hand,  equation  (B5)  can  be  obtained  from  equation  (98), 
if  m(x)  =  m,  n(x)  =  n,  and  k' (x)  =  are  substituted  and  the 

O 

b  term  is  disregarded.  Thus,  in  the  limiting  case  considered,  Ferri’s 
and  the  present  results  are  identical. 
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APPENDIX  C 

AIRFOIL  GEOMETRY 


In  order  to  show  that  the  considerations  in  the  text  pertain  to 
airfoils  at  both  zero  and  finite  angles  of  attack,  the  geometry  of  the 
airfoil  will  be  considered  in  detail. 

If  the  chord  of  the  airfoil  is  at  zero  angle  of  attack,  the  local 
slope  is  given  by 


=  ek'(x)  Z  e©  (Cl) 

dx 


where  y  and  k(x)  are  positive  on  the  upper,  and  negative  on  the 
lower  surface  (see  fig.  1(a)).  If  the  airfoil  is  symmetrical,  then: 


Figure  l(b)  may  illustrate  how  these  relations  have  to  be  modified 
in  case  the  chord  of  the  airfoil  is  at  a  finite  angle  e6  relative  to 
the  flight  direction  x. 


Two  coordinate  systems  may  be  used,  one  whose  x-axis  coincides  with 
the  flight  direction  and  another  one  whose  x-axis  coincides  with  the 
chord.  Then  the  slope  of  the  profile  above  the  flight  direction  x  is: 


d y 

dx 


ek'  (x) 


=  tan  e(@  -  9) 


e(6  -  0) 


(C3) 
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where  tan  is  the  slope  of  the  profile  above  the  chord.  If  the 
profile  function  above  (or  below)  the  chord,  y  =  ek(x),  is  introduced, 
then  ’ 


ee  ~  tan  eg  =  k*  (x) 


and  equation  (C3)  can  be  rewritten: 


(c4) 


Note  that  y  >  0  on  the  upper  surface  and  y  <  0  on  the  lower  surface. 
Then,  at  P  _in  figure  l(b),  e@  is_positive  because  y  increases  with 
increasing  x.  Fixing  the  sign  of  9  analogously  to  that  of  6,  ~Q  in 
the  sketch  is  positive  and  larger  than  9,  so  that  in  equation  (C3) 

'  y 

comes  out  negative,  as  it  should  according  to  the  figure. 

Now,  the  barred  and  unbarred  coordinate  systems  are  related  (neg¬ 
lecting  second  orders  in  e)  as ’follows: 


x  ~  x  -  y0 

y  ~  y  +  x0 


Then 


ek1  (x)  =  ek’  [x(l  +  e2S2)  -  9  k(x)e^J 


ek’ (x) 


Therefore,  equation  (c4)  may  be  written: 


dy 

dx 


ek’ (x) 


(C5) 


\ 


e|k  ’(x)  -  9] 
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Figure  3*-  Total  drag  coefficient  Cp  as  function  of  Mach  number 

for  diamond  profile,  e  =  0.05* 


for  diamond  profile,  g  =  0.05. 
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Figure  7-_  Total  pitching-moment  coefficient  Cj^  as  function  of  Mach 
number  M0  for  diamond  profile.  e  =  0.05- 
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